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ABSTRACT 

We propose that it should be possible to use the CMB to discriminate between 
dark energy models with different equations of state, including distinguishing a 
cosmological constant from many models of quintessence. The separation of peaks 
in the CMB anisotropics can be parametrised by three quantities: the amount of 
quintessence today, the amount at last scattering, and the averaged equation of 
state of quintessence. In particular, we show that the CMB peaks can be used 
to measure the amount of dark energy present before last scattering. 

Subject headings: cosmic microwave background — cosmology: theory 



1. Introduction 



The idea of quintessence was born (Wetterich 1988) from an attempt to understand the 
vanishing of the cosmological constant. It was proposed that the cosmological evolution of a 
scalar field may naturally lead to an observable, homogeneous dark energy component today. 
This contrasts with the extreme fine-tuning needed in order for a cosmological constant to 
become significant just at recent times. If quintessence constitutes a major part of the energy 
density of the Universe today, say > 0.5, structure formation tells us that this cannot 
always have been so in the past (Peebles & Ratra 1988; Ratra & Peebles 1988; Ferreira & 
Joyce 1997, 1998). Combining the phenomenology of a large quintessence component with 
the quest for naturalness (Hebecker & Wetterich 2000) leads to cosmologies with an equation 
of state for quintessence changing in time, compatible with a universe accelerating today. 
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In several aspects of phenomenology the models with a dynamical dark matter component 
resemble the cosmology with a cosmological constant (Huey et al. 1999). It is therefore 
crucial to find possible observations which allow us to discriminate between the dynamical 
quintessence models and a constant dark energy theory (i.e. a cosmological constant). The 
detailed structure of the anisotropies in the cosmic microwave background radiation (CMB) 
depends upon two epochs in cosmology: around the emission of the radiation (last scattering) 
and today. The CMB may therefore serve as a test to distinguish models where quintessence 
played a role at the time of last scattering from those where it was insignificant at this 
epoch. It may also reveal details of the equation of state of quintessence (characterised by 
w = p^/p^) in the present epoch. 

The calculation of CMB spectra is, in general, an elaborate task (Seljak & Zaldarriaga 
1996; Hu & Sugiyama 1995). However, the location of the peaks and, for our purpose, the 
spacing between the peaks can be estimated with much less detailed knowledge if adiabatic 
initial conditions and a flat universe are assumed. The oscillations of the primeval plasma 
before decoupling lead to pronounced peaks in the dependence of the averaged anisotropies 
on the length scale. When projected onto the sky today, the spacing between the peaks at 
different angular momentum / depends, in addition, on the geometry of the universe at later 
time. It is given, to a good approximation, by the simple formula (Hu & Sugiyama 1995; 
Hu et al. 1997) 

a r T Q -n s T -r is 

Al = 71 = 7T— . (1) 

S C s Tis 

Here Tq and T\ s are the conformal time today and at last scattering (which are equal to 
the particle horizons) and r = J dt a _1 (t), with cosmological scale factor a. The sound 
horizon at last scattering s is related to T\ s by s = c s T\ s , where the average sound speed 
before last scattering c s = r^ 1 JJ" 15 dr c s obeys c~ 2 = 3 + (9/4)(p b (t)/p r (t)), with p b /p r the 
ratio of baryon to photon energy density. We note a direct dependence of Al on the present 
geometry through r as well as an indirect one through the dependence of T; s on the amount 
of dark energy today (see Equation (7)). 

The location of the m— th peak can be approximated by (Hu et al. 2000) 

l m = Al (m - <j>) , (2) 

where the phase-shift <fi is typically less than 0.4 and is determined predominantly by re- 
combination physics. By taking the ratio of two peak locations (say Z1//2), the factor Al 
and with it the dependence on post-recombination physics drops out and we are in principle 
able to probe pre-recombination dark energy directly. If the other cosmological parameters 
were known the dependence of on the amount of dark energy at last scattering could pro- 
vide a direct test of this aspect of quintessence models. Unfortunately, <fi also depends on 
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other cosmological parameters including the baryon density and spectral index and there is 
no known analytic formula for <fi (and in fact does have some m— dependence). We first 
concentrate on the peak spacing Al for which an analytic formula can be given. 

The equation of state of a hypothetical dark energy component influences the expansion 
rate of the Universe and thus the locations of the CMB peaks (Ferreira & Joyce 1997, 1998; 
Coble et al. 1997; Caldwell et al. 1998; Amendola 2000). In particular the horizons at 
last scattering and today are modified, leaving an imprint in the spacing of the peaks. The 
influence of dark energy on the present horizon and therefore on the CMB has been discussed 
in (Brax et al. 2000). A likelihood analysis on combined CMB, large scale structure and 
supernovae data (Efstathiou 1999; Bond et al. 2000) can also give limits on the equation of 
state. Several of these analysis concentrate on models where the dark energy component is 
negligible at last scattering. In contrast, we are interested particularly in getting information 
about dark energy in early cosmology. Therefore, the amount of dark energy at last scattering 
is an important parameter in our investigation. 




Fig. 1. — The CMB Spectrum for A-CDM (model C), leaping kinetic term (model A) and 
inverse power law (model B) quintessence universes with = 0.6. The data points from 
the Boomerang (de Bernadis et al. 2000) and Maxima (Hanany et al. 2000) experiments are 
shown for reference. 
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We present here a quantitative discussion of the mechanisms which determine the 
spreading of the peaks. A simple analytic formula permits us to relate Al directly to three 
characteristic quantities for the history of quintessence, namely the fraction of dark energy 
today, Qq, the averaged ratio between dark pressure and dark energy, Wq = (p^/p^)o, and 
the averaged quintessence fraction before last scattering, Clf s (for details of the averaging 
see below). We compare our estimate with an explicit numerical solution of the relevant 
cosmological equations using CMB-FAST (Seljak & Zaldarriaga 1996). For a given model of 
quintessence the computation of the relevant parameters SIq , Wo and Q ls requires the solution 
of the background equations. Our main conclusion is that future high-precision measure- 
ments of the location of the CMB-peaks can discriminate between different models of dark 
energy if some of the cosmological parameters are fixed by independent observations. It 
should be noted here that a likelihood analysis of the kind performed in (Bond et al. 2000), 
where w is assumed to be constant throughout the history of the Universe, would not be 
able to extract this information as it does not allow £lf s to vary. We point out that for 
time- varying w there is no direct connection between the parameters Wo and Q ls , i.e. a 
substantial Q ls (say 0.1) can coexist with rather large negative Wq. We perform therefore a 
three parameter analysis of quintessence models and our work goes beyond the investigation 
for constant w in (Huey et al. 1999). 



2. CMB Peaks in Quintessence Models 



We wish first to illustrate the impact of different dark energy models on the fluctuation 
spectrum of the CMB by comparing three examples. The first corresponds to a 'leaping 
kinetic term quintessence' (Hebecker & Wetterich 2000) (A), the second to 'inverse power- 
law quintessence' (Peebles & Ratra 1988; Ratra & Peebles 1988) (B) and the third to a 
cosmological constant (C). The three examples, whose parameters are chosen such that = 
0.6 for each, give similar predictions for many aspects of cosmological observation (we assume 
everywhere a flat universe Jltotai — !)• Details of the models can be found below in Section 
5. We solve the cosmology using CMB-FAST for a flat initial spectrum with parameters 
specified in Table 3. Figure 1 clearly demonstrates that the fluctuation spectra of the three 
models are distinguishable by future high-precision measurements. This can be traced back 
to different values of Qf s and w , namely Qf s = (0.13, 0, 0) and w = (—0.45, —0.37, —1) for 
models (A,B,C). These quantities enter a simple analytic formula (derived below in Section 
3) for the spacing between the peaks 



Al = nc: 




(3) 
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with 

and today's radiation component Q T = 9.89 x 10 -5 , a^ 1 = 1100 and c s = 0.52. In Table 
1 we evaluate Equation (3) for quintessence models with various parameters, together with 
the locations l 1: l 2 of the first two peaks computed by CMB-FAST. The last entry contains 
the peak spacing averaged over 6 peaks for the numerical solution. This demonstrates that 
an accurate measurement of the peak spacing Al is a powerful tool for the discrimination 
between different dark matter models! 



3. Analytic Estimate of Peak Spacing 

We derive next the formula (3). Our first task is to estimate the sound horizon at 
decoupling. We assume that the fraction of quintessential energy Q^(t) does not change 
rapidly for a considerable period before decoupling and define an effective average Q ls = 
T is 1 Jo^ 8 ^( r )dT. We note that this average is dominated for r near T\ s whereas very early 
cosmology is irrelevant. Approximating £7^ by the constant average for the period around 
last scattering, the Friedmann equation for a flat universe reads 



3M?tf 2 (t)(l - nj) = p m (t) + p r (t) = p™a(t)- 3 + pla(ty 



(5) 



Here Mp = (SttGnY" 1 ^ 2 ^ is the reduced Planck mass, H(t) is the Hubble parameter and p™ 
and Pq are the matter and relativistic (photons and 3 species of neutrinos) energy densities 
today. 



Today, neglecting radiation, we have ?>M"pHl{l — VLq] 
(5) to obtain 



p™, which we insert in Equation 



da V 
drj 



tfoa-^itr 1 (1 - nj)a(r) + nj 



(6) 



where we have changed from coordinate time t to conformal time r. Separating the variables 
and integrating gives 




(7) 



which is well known for vanishing Q h . For fixed H Q , Qq, Q t and a\ s (see Table 3 for the 
values used in this paper), we see that T\ s = (1 — f2^)( 1//2 \ where r™ is the last scattering 



- 6- 



horizon for a A-CDM universe (which we treat here to be just a special realisation of dark 
energy with w — — 1). To estimate the sound horizon, we also need c s , which may be obtained 
numerically and in our model universe is 0.52. 

Turning to the horizon today, we mimic the steps of above, this time assuming some 
equation of state p^(t) = w(t)p^(t) for quintessence. 

We define an averaged value Wq by 

/TO / PTo \ -1 

n*(r)iu(r)dr x fy ^(r)drj . (8) 

It is fi^-weighted, reflecting the fact that the equation of state of the dark energy component 
is more significant if the dark energy constitutes a higher proportion of the total energy of 
the Universe (see Figure 2). 

In the limiting case that the equation of state did not change during the recent history of 
the Universe, the average is of course equal to w today. Nevertheless, the difference between 
the average Wo and today's value wq can be substantial for certain models, as can be seen 
from Table 2. 

Integrating the cosmological equation with constant w 

(l?) 2 = H ° { (1 ~ "° ~ QT ° Ht) + ^ ail ~ 3m) + Q o} ' ( ) 

gives 

r = 2H \l-nt)^F(niw ), (10) 

with F given by Equation (4). Substituting Equations (7) and (10) into Equation (1), we 
obtain the final result (3). 

The integral F of Equation (4) can be solved analytically for special values of w , e.g. 

F(a*,wo = 0) = y/l-Slt (l - v^f) + 0(n&). (11) 

Since the integral (4) is dominated by a close to one (typically w < 0) only the present 
epoch matters, consistent with the averaging procedure (8). From this we regain on inserting 
in Equation (10) the trivial result that the age of the Universe is the same for a cold dark 
matter and a pressureless dark energy universe. We plot F(Qq,w ) for various values of Qq 
in Figure 3. 

For < 0.6, Equation (3) to good approximation (better than one percent) can be 
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Fig. 2. — Equation of state w(t), w(t) = VI^(t)w(t)tq/ J q t ° Q^(t')(\t' and averaged equation 
of state Wq for the leaping kinetic term model with ilf s = 0.13 and £lf s = 0.22. 



written 



A/ = tic- 1 



F(njw ) 
a is (l- ilt 



1/2 



fi r n 



a\ s ( 1 - fir 



2a ls 1 - fi 



(12) 



The precision of our analytic estimate for Al can be inferred from Table 1. Similarly, 
we show in Table 2 the accuracy of the estimates of r ls (7) and r (10) by comparison 
with the numerical solution. This demonstrates that our averaging prescriptions are indeed 
meaningful. We conclude that the influence of a wide class of different quintessence models 
(beyond the ones discussed here explicitly) on the spreading of the CMB-peaks can be 
characterised by the three quantities fig, fi ls and wq. 



4. Ratios of peak locations 



An alternative to the spacing between the peaks is the ratio of any two peak (or indeed 
trough) locations. After last scattering the CMB anisotropics simply scale according to the 
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Fig. 3. — F(Qq,Wo) as a function of Wo of the dark energy component, for Oq between 
0.2 and 0.7. Between the limiting cases of Wo = —1 (cosmological constant) and Wo = 
(corresponding to pressureless dust), the age of the Universe varies considerably. 



geometry of the Universe - taking the ratio of two peak locations factors out this scaling 
and leaves a quantity which is sensitive only to pre-last-scattering physics. As can be seen 
in Table 1, (spatially-flat) models with negligible ilf s all have l 2 /h ~ 2.41 for the parameters 
given in Table 3. The dependence of this ratio on the other cosmological parameters can 
be computed numerically (Doran & Lilley 2001). If the other parameters can be fixed 
by independent observations, the ratios of peak locations are fixed uniquely for models 
with vanishing ilf s . A deviation from the predicted value would be a hint of time- varying 
quintessence. It may also be possible to make a direct measurement of fi^ from ratios of 
successive peak locations. 
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5. Specific Quintessence Models 

Different models of quintessence may be characterised by the potential V(<f>) and the 
kinetic term of the scalar 'cosmon'-field 

m = l -{d,<pfk 2 (<p) + v{ ( i ) ). (13) 

For practical purposes, a variable transformation allows us to work either with a standard 
kinetic term k(4>) = 1 or a standard potential, i.e. V(<f>) = Mp exp(— <p/Mp). A cosmological 
constant corresponds to the limit V((f>) = A, k((f>) = 0. It is also mimicked by k((f>) — > oo. 
We consider four types of model. 

A. A 'leaping kinetic term' model (Hebecker & Wetterich 2000), with 

V{<j>) = V(<f>) = M 4 p exp(-<f)/Mp) (14) 

and kinetic term 

k (0) = k min + tanh [(0 - fa) /Mp] + 1. (15) 

We have taken k m \ n = 0.05, 0.1, 0.2 and 0.26 and <p\ is adjusted to ~ 277 in order to 
obtain = 0.6. The value of flf s is determined by these parameters. 

B. An inverse power law potential (Peebles & Ratra 1988; Ratra & Peebles 1988), with 

k(4>) = 1 and 

Vtf) = A<j>- a , (16) 

We have chosen a = 6, 22 and 40, and A adjusted such that = 0.6. Once again, Q ls 
follows. 

C. A cosmological constant tuned such that Qq = Qq = 0.6. 

D. The original exponential potential (Wetterich 1988, 1995), with k{4>) = 1 and 

V{(j>) = M|exp(- v / 2a0/Mp), (17) 



where a = 




For the models (A) and (D), quintessence is not negligible at last scattering. The pure 
exponential potential requires Qq < 0.2 for consistency with nucleosynthesis and structure 
formation. It does not lead to a presently accelerating universe. We quote results for 
nf; = o.6 for comparison with other models and in order to demonstrate that a measurement 
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of Al can serve as a constraint for this type of models, independently of other arguments. 
The inverse power law models (B) are compatible with a universe accelerating today only 
if fi^ is negligible. Again, our parameter list includes cases which are not favoured by 
phenomenology. As an illustration we quote in Table 1 the value of <7g, which should typically 
range between 0.6 and 1.1 for the models considered. For example, the exponential potential 
model with large is clearly ruled out by its tiny value of erg 1 . The main interest for 
listing also phenomenologically disfavored models arises from the question to what extent 
the location of the peaks can give independent constraints. From the point of view of 
naturalness, only the models (A) and (D) do not involve tiny parameters or small mass 
scales. 

The horizons and Al for the models considered are shown in Tables 2 and I. We note 
that the estimate and the exact numerical calculation are in very good agreement. A different 
choice of a\ s , say c^" 1 = 1150, would have affected the outcome on the low-percent level. Also, 
the average spacing obtained from CMB-FAST varies slightly (at most 2%) when averaging 
over 4, 5 or 6 peaks. For a fixed value of the equation of state, Wq = —0.7, we plot the peak 
spacing as a function of fig and Q,f s in Figure 4. 

For fixed w and fig, we see from Equation (3) that Al oc (1 — fi^)^ 1 / 2 ). Hence, when 
combining bounds on Qq and Wo from the structure of the Universe, supernovae redshifts 
and other sources with CMB data, the amount of dark energy in a redshift range of z ~ 10 5 
to last scattering z ~ 1100 may be determined. 

From Figure 2, we see that the averaged equation of state of the quintessence field for the 
present epoch is, in principle, a very influential quantity in determining the spreading of the 
peaks. Since combined large scale structure, supernovae and CMB analysis in (Bond et al. 
2000) suggest Wq < — 0.7, the difference between a cosmological constant and quintessence 
may be hard to spot if £lf s is negligible. However, even with the data currently available, the 
first peak is determined to be at I = 212 ± 7 (Bond et al. 2000). Once the third and fourth 
peak have been measured, the measurement of the spacing between the peaks becomes an 
averaging process with high precision. We can then hope to distinguish between different 
scenarios. 



1 Of course cr$ itself also depends on other cosmological parameters and so it alone cannot be used to 
determine fl ls . 



0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 

Is 



Fig. 4. — Contours of equal peak spacing A/ as a function of and f2 ls .The average 
equation of state is kept fixed, Wq = —0.7. Increasing nf s leads to a pronounced stretching 
of the spacing. 



6. Conclusions 



The influence of quintessence on the spacing between the CMB peaks is determined 
by three quantities: SIq , Qf s , and w . When the location of the third peak is accurately 
measured, we can hope to be able to discriminate between a pure cosmological constant and 
a form of dark energy that has a non-trivial equation of state - possibly, and most likely, 
changing in time. The peak ratios will help determining which in principle can also be 
extracted from A/, if and Wq are measured by independent observations. With Q ls fixed, 
the peak spacing can be used to constrain fl^ and wo- This can permit consistency checks for 
the quintessence scenario. Together with bounds on fl^ for the period of structure formation 
(5 < z < 10 4 ) and the bound &bbn < 0-2 from big bang nucleosynthesis (Wetterich 1988, 
1995; Birkel & Sarkar 1997) (z ~ 10 9 ) we will post a few milestones in our attempt to trace 
the cosmological history of quintessence. 
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Table 1. Location and spacing of the CMB peaks for several models a 



nl w h l 2 l 2 /h A/ cstim - Al 



Leaping kinetic term (A), fig = 0.6 



8.4 x 10~ 3 


-0.76 215 518 2.41 


292 


291 


0.86 


0.03 


-0.69 214 520 2.43 


294 


293 


0.78 


0.13 


-0.45 211 523 2.48 


299 


300 


0.47 


0.22 


-0.32 207 524 2.53 


302 


307 


0.29 


Inverse power law potential (B), Qq = 0.6 


8.4 x 10~ 8 


-0.37 199 480 2.41 


271 


269 


0.61 


9.9 x KT 2 


-0.13 178 443 2.49 


252 


252 


0.18 


0.22 


-8.1 x 10~ 2 172 444 2.58 


257 


257 


0.09 




Pure exponential potential, Og 


= 0.6 






0.70 


7 x 10~ 3 190 573 3.02 


368 


377 


0.01 




Pure exponential potential, Qq 


= 0.2 






0.22 


4.7 x 10~ 3 194 490 2.53 


282 


281 


0.38 




Cosmological constant (C), Qq 


= 0.6 









-1 219 527 2.41 


296 


295 


0.97 




Cold Dark Matter - no dark energy, fig = 










205 496 2.42 


269 


268 


1.49 



a The analytic estimate of Al stems from Equation (3). The position 
of the peaks l\ and l 2 and A/ num - are calculated with a modified CMB- 
FAST code. 
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Table 2. Horizons in Mpc at last scattering and today for various kinds of quintessence. a 



Ql Wq w r cstim - r num - Ar r, e s stim - Ar is 



Leaping kinetic term (A), Oq = 0.6 



8.4 x 10- 3 


-0.79 -0.76 13073 13060 


0.1% 


266 


0.3% 


0.03 


-0.79 -0.69 12971 13000 


0.2% 


263 


0.3% 


0.13 


_n 78 _n 45 12470 12590 


1 0% 

J. . \j / u 


248 


2% 


0.22 


-0.75 -0.32 12012 12175 


1.3% 


236 


0.0% 




Inverse power law potential (B), Qq 


= 0.6 






8.4 x 10~ 8 


-0.32 -0.37 12205 12147 


0.5% 


267 


0.0% 


9.9 x 10~ 2 


-0.16 -0.13 10774 10798 


0.2% 


253 


0.2% 


0.22 


-0.1 -8.1 x 10~ 2 10241 10273 


0.3% 


236 


0.2% 




Pure exponential potential, £Iq = 


0.6 






0.70 


0.00 7 x 10~ 3 9014 9049 


0.4% 


146 


2.3% 




Pure exponential potential, Qq = 


0.2 






0.22 


5 x 10~ 5 4.7 x 10~ 3 9107 9120 


0.1% 


191 


0.3% 




Cosmological constant (C), = 


0.6 









-1 -1 13330 13325 


0.0% 


267 


0.0% 




Cold Dark Matter - no dark energy, 


nt = 









9133 9133 


0.0% 


201 


0.5% 



a We compare analytical estimates with numerical solutions. We also list the 
equation of state today wq (relevant for supernovae observations) and compare 
it with the averaged equation of state Wo- The cosmological parameters used 
can be read off Table 3 . 



-16- 



Table 3. Symbols, their meanings and numerical values used in this paper. 



Symbol 


Meaning 


Value 


M P 


reduced Planck mass Mp 2 = 8irG 




a(r) 


scale factor, normalised to unity today 




a\ s 


scale factor at last scattering 


noir 1 


h 


Hubble parameter today H = 100/iokm s _1 Mpc _1 


0.65 


n 


relativistic Q today 


9.89 x 10~ 5 




baryon Q today 


0.05 


C s 


r-averaged sound speed until last scattering 


0.52 


n 


spectral index of initial perturbations 


1 



